order=4
period=3

fom_pos=fom_pos.model
fom_neg=fom_neg.model
n=1000
length=90
maxlen=300

train_pos=pos.train
test_pos=pos.test
train_neg=neg.train
test_neg=neg.test

trained_fom_pos=trained.fom_pos.model
trained_vom_pos=trained.vom_pos.model
trained_fom_neg=trained.fom_neg.model

likelihood_pos_pos=test_pos.pos.ll
likelihood_pos_neg=test_pos.neg.ll
likelihood_neg_pos=test_neg.pos.ll
likelihood_neg_neg=test_neg.neg.ll

test_acc=test.acc
test_stats=test.stats
test_result=test.result

cat /dev/null > $test_result
etas="0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1"
for eta1 in $etas
do 
eta=0$(echo "$eta1 / $period" | bc -l)
./vom --rand_fom --order=4 --period=$period --output=$fom_pos
./vom --rand_fom --order=0 --period=$period --output=$fom_neg

./vom --simulate --model=$fom_pos --output=$train_pos --n=$n --length=$length --period=$period --maxlen=$maxlen
./vom --simulate --model=$fom_pos --output=$test_pos --n=$n --length=$length --period=$period --maxlen=$maxlen
./vom --simulate --model=$fom_neg --output=$train_neg --n=$n --length=$length --period=$period --maxlen=$maxlen
./vom --simulate --model=$fom_neg --output=$test_neg --n=$n --length=$length --period=$period --maxlen=$maxlen


./vom --fom --input=$train_pos --order=$order --output=$trained_fom_pos --period=$period
./vom --grow --input=$train_pos --model=$trained_fom_pos --output=$trained_vom_pos --eta=$eta --prune --period=$period

./vom --fom --input=$train_neg --order=0 --output=$trained_fom_neg --period=$period

./vom --likelihood --input=$test_pos --model=$trained_vom_pos --output=$likelihood_pos_pos --period=$period
./vom --likelihood --input=$test_pos --model=$trained_fom_neg --output=$likelihood_pos_neg --period=$period
./vom --likelihood --input=$test_neg --model=$trained_vom_pos --output=$likelihood_neg_pos --period=$period
./vom --likelihood --input=$test_neg --model=$trained_fom_neg --output=$likelihood_neg_neg --period=$period

./vom --accuracy --negneg=$likelihood_neg_neg --negpos=$likelihood_neg_pos --pospos=$likelihood_pos_pos --posneg=$likelihood_pos_neg --output=$test_acc.$eta
./vom --depth --model=$trained_vom_pos --period=$period --output=$test_stats.$eta

echo $eta1 `head -n 1 $test_acc.$eta` `head -n 1 $test_stats.$eta` >> $test_result

done
